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The Josephson-like interband couplings in multi-band superconductivity exhibit degenerate energy minima, 
which support states with kinks in phase of superconductivity. When the interband couplings in systems of 
three or more components are frustrated, the time-reversal symmetry (TRS) can be broken, which generates 
another type of phase kink between the two time-reversal-symmetry breaking (TRSB) pair states. In this work, 
we focus on these novel states of phase kinks, and investigate their stability, similarity, differences and physical 
consequences. Main results are summarized as follows: (1) We find a new type of phase slip when the kink 
becomes unstable. (2) In the kink region, TRS is broken and spontaneous magnetic fields are induced. (3) In 
superconductors with TRSB, composite topological excitations associated with variations of both superconduc- 
tivity phase and amplitude can be created by local perturbations, or due to proximity effect between normal 
metals. 
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Introduction - It has been known for long time that supercon- 
ductors with different pairing symmetries in contact with one 
another can form stable domain structures HI |2j. Properties of 
domain walls are governed by the pairing symmetries in the 
domains, thus these heterogeneous systems become vital in 
understanding the pair symmetry. Meanwhile, there are grow- 
ing evidences that superconductors may break discrete sym- 
metries in addition to the U(l) (local) gauge symmetry whose 
loss defines superconductivity [3, 4j. Examples include TRSB 
in some unconventional superconductorsfS] [6), which results 
in unusual phenomena such as the appearance of magnetic 
flux when the superconductivity is perturbed by nonmagnetic 
impurities|7|. These superconductors can also form stable do- 
main walls between domains of distinct symmetry-breaking 
states. 

The discovery of MgB2|E] and iron-pnictide 
superconductors (9) has opened intense and exciting dis- 
cussions of multi-band superconductivity in condensed 
matter physics. In these systems, superconductivity in one 
band is coupled through interband Josephson coupling to that 
in another band y, / A,A / cos(0; - <f>j), with tp, and A,- being 
the superconductivity phase and amplitude in the i-th band 
respectively. This gives rise a collective oscillation of the 
superconductivity phases, known as the Leggett mode| 10 1. 

It is interesting to observe that the interband coupling has 
degenerate energy minima <p, - tpj = 2nn for < 0, which 
supports various topological excitations in the form of phase 
kinks belonging to the homotopy class 7To{S°), whereas the 
well known vortex solution in type II superconductors belongs 
to the homotopy class ni(S '). The existence of the kink so- 
lution was first discussed by Tanaka| 11 1 for two-component 
superconductors, and later it was discussed that phase kinks 
can be excited in nonequilibrium processes such as current 
iniection |[T2l . The phase kinks have been observed experi- 
mentally in layered aluminium mesocopic rings with two or- 
der parameters [ 1 3 1 . 

In the presence of frustrated interband couplings in super- 



conductors with three or more components, the system may 
break the TRS lfl4l4T8ll . In the TRSB state, $ + e ie W* for 
any phase 8 with 4* = (¥1,^2, — ,*Pn) a vector of the com- 
plex order parameters. A phase kink may appear between two 
degenerate states 4* and T*. One thus sees that multi-band su- 
perconductors support two types of kink solutions of different 
origins. In a recent paper by Garaud et. al. [19], composite 
topological excitations associated with phase kink and vortex 
in superconductors with TRSB have been found numerically. 
The stability of these phase kinks however still remains to be 
investigated. 

In the present work, we investigate the stability, similar- 
ity, differences and physical consequences of these two types 
of kink solutions. In superconductors with TRS in bulk, the 
phase kink breaks TRS at the domain wall, and induces local 
magnetic flux. Upon elevation of temperatures, kinks become 
unstable as a consequence of increasing coherence length. At 
the instability, a phase slippage occurs accompanying a volt- 
age pulse. Contrarily, kinks between TRSB pair states remain 
stable even with the increasing coherence length. Moreover, 
in superconductors with TRSB, various types of composite 
topological excitations associated with the variation of super- 
conductivity phase and amplitude can be created by perturb- 
ing the superconductivity locally, such as heating and/or non- 
magnetic impurities, and at interface to a normal metal. 

Kink solutions - We start from the standard multi-band 
Ginzburg-Landau (GL) theory with Josephson-like interband 
couplings ll20l [21], which is adequate for discussions on 
physics addressed here 
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where symbols are conventionally defined El . Throughout 
the paper, we use the units % = 2e = c = 1. ju for I + j is the 
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FIG. 1. (color online). Domain structure in multi-band supercon- 
ductors (left), phase kink in a two-band superconductor (middle), and 
that in a three-band superconductor with each domain corresponds to 
distinct TRSB pair states (right). 

interband coupling, which can be either repulsive or attrac- 
tive depending on the strengths of the Coulomb and electron- 
phonon interactions. The interband repulsion may cause frus- 
tration of the superconductivity in different bands and results 
in TRSB lfl4m8l . For MgB 2 , the interband coupling is com- 
monly accepted as attractive 712 < 0, while for iron-pnictide 
superconductors, there are growing evidences that some of y^ 
are positive and the system favors s+ pair symmetry. Il23l l24l 
The kinetic energy in Eq. ([T]i can be rewritten as 

4^ |(V + iA^jf, where ^ = -^1 ]2mj\aJ\ and H cj are co- 
herence length and thermodynamic critical field in respective 
single-band condensates (yy = 0). When width of the kink 
Ak (derived below) is much larger than gj, Ak » £/, the sup- 
pression of the amplitude of the order parameters by the phase 
kink is weak, and the order parameter is approximately con- 
stant in space. In this case, we can concentrate on the phase 
variables of the order parameters. 

First we consider phase kink between domains with TRS in 
one dimension, where we can take the gauge A = 0. The min- 
imal model for this domain structure is of two bands. Since 
the sign of yn can be gauged away in this case, we consider 
j\2 < without loss of generality. We also assume an identi- 
cal amplitude of order parameter A, = A for simplicity. The 
variation of the phase difference <pi2 = <P\ - fyi is described by 
the sine-Gordon equation lfTTl 

d 2 K 4>u + 2j\i(m x + m 2 ) sin 4> n = 0, (2) 

and d x (f>i = -m\d x (p2l{m\ + mi)- The width of the kink is 
Ak — 1 / -y/-2yi2(»ii + f«2), which is temperature independent. 
The condition that A^ » £ ; then becomes \yn\ <k or,-. A typical 
phase kink is shown in Fig. [T] (middle). The TRS is broken 
at the domain wall while reserved in the domains. There are 
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FIG. 2. (color online). Numerical results of the magnetic field 
distribution: (a) a circular domain wall in a two-band superconduc- 
tor; (b) a circular domain wall between two TRSB pair states in a 
three-band superconductor. For (a), (Xj = -20, 712 = — 1, fij = 1, 
mi = 1 and m 2 = 3 in the numerical calculations; for (b), aj = 0, 
j3j = irij = pj = 1, 712 = 1, 713 = 1.2 and 723 = 1.5. 

finite phase differences between the right and left domains in 
both components. 

Now we consider the kink solution of a superconductor 
of three or more components with frustrated interband cou- 
plings, where TRS is broken in bulk. As a minimal model, 
we treat a superconductor with three equivalent bands aj = a, 
jij — y > 0, and /«, = m. The two degenerate ground states 
W = A(l,e i2 " /3 ,e i4 ' l/3 ) and = A(l,e- i2K/ \e- i4 " /3 ) as a con- 
sequence of TRSB are displayed in Fig[T](right). For constant 
amplitudes of order parameters at y <K a, the phase kink is 
described by d x <p\ = 0, 3 v (0i2 + ^13) = and 

1 , 

- — d-(f>n + sin 0i2 + sin (2<j> n ) = 0. (3) 
Zmy 

The potential associated with Eq. d3J V p = cos0i 2 + 
cos(20i2)/2 has many degenerate minima (f>i2, m = ±2n/3 + 
2nn. One can construct kink solution between any pair of 
the energy minima. Their stability and magnetic response are 
qualitatively the same. To be specific, we only consider the 
following kink solution, which can be found analytically us- 
ing the Bogomolny inequality 1 25 1 
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and the associated energy is 

E k = 3 V^(3V3-;r). (5) 

In one dimension (ID), there is no supercurrent in the do- 
main wall due to the current conservation d x J s = 0. In higher 
dimensions, supercurrent and the associated magnetic field are 
induced at the domain wall as a result of TRSB. We consider 
a closed domain wall described by either Eq. |2]i or Eq. (|4]) 
in a 2D superconductor. To investigate the dynamic evolution 



3 



of the domain wall, we solve the time-dependent GL equa- 
tion(TDGL) numerically! 26 1 
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with Dj the diffusion constant, cr the normal conductivity, and 
O the electric potential. 

In simulations, we prepare a closed domain wall with 
square or rectangular shapes as initial conditions. In order 
to minimize its energy, the domain wall organizes itself into 
a circular shape irrespective to its initial shape during the 
time evolution in simulations. Magnetic fields appear spon- 
taneously at the domain wall with alternating directions, as 
shown in Fig. [2] (a) and (b). As revealed by numerical sim- 
ulations, for phase kinks in superconductors with TRS [see 
Eq. Q], the induced magnetic field changes polarization in 
both radial and azimuthal directions as shown in Fig. [21(a). 
While for kinks between TRSB pair states [see Eq. d4b]7 the 
magnetic field changes polarization only in the azimuthal di- 
rection as shown in Fig. |2](b). One may treat the domain wall 
at the left semicircle as a phase kink, then the domain wall at 
the right semicircle is an anti-kink. They attract each other, 
which causes the whole circular domain wall collapsing, and 
renders a uniform state. Since the attraction between two do- 
main walls becomes exponentially weak at a large separation, 
the life time of the domain walls increases with the size of the 
domain enclosed. This allows for possible experimental de- 
tections on the induced magnetic flux after quenching when 
domain walls are excited by chance. 

Stability of the Kink Solution- We proceed to investigate the 
stability of the kink solution in Eq. (|2| taking into account 
the suppression of amplitude of order parameter by the phase 
kink. The magnitude of the suppression depends on the ratio 
of the kink width to the coherence length £ as briefly men- 
tioned above. As the coherence length increases when tem- 
perature is elevated while the width of kink remains almost 
unchanged, the superconductivity in the domain wall will be 
greatly depleted. At a threshold value, the phase kink loses 
its stability, and system evolves into a uniform state. There is 
a voltage pulse associated with varying magnetic field across 
the domain wall which is experimentally detectable. This pro- 
cess is a new type of phase slip, different from that in single- 
band superconductors carrying supercurrent close to the criti- 
cal one, with the latter one caused by fluctuations. [22 1 

We explicitly consider a superconductor of two identical 
bands with a phase kink localized at the center of a super- 
conducting wire. We solve numerically the TDGL equations 
and derive the stable configuration of the superconductivity 
phase as temperature (namely a) varies. When temperature 
increases, the amplitude of the superconductivity at the do- 
main wall decreases as depicted in Fig. |3ja). At a threshold 



a for given value of j\i [symbols in Fig.|3jb)], the phase kink 
becomes unstable and the system evolves into a uniform state, 
during which a voltage pulse appears. Therefore, the phase 
kinks in superconductors with TRS are stable only for weak 
interband couplings. 

A superconducting wire with a phase kink can be alterna- 
tively considered as a Josephson junction since the supercon- 
ductivity is suppressed at the domain wall. In the ground state, 
the phase difference between two domains is finite, thus it is 
a realization of 0-junction[27 1, or 7r-junction[28] if the two 
bands are identical. When current is injected into the wire, 
the phase kink is deformed due to the phase gradient created 
by the injected current. At a threshold current, the phase kink 
becomes unstable, and the system evolves into a uniform su- 
perconducting state, during which a voltage pulse appears, 
similar to the case with increasing temperature. The thresh- 
old current is still much smaller than the depairing current of 
the uniform state. One may regard the threshold current as a 
critical current for the present Josephson junction. 

We perform numerical calculations on the critical current of 
a superconducting wire with a phase kink, introducing super- 
current into the system by twisting the phases at the two edges 
of the wire far away from the phase kink. The kink structure 
is deformed into the shape depicted in Fig. [3f c) by the current 
injection. The critical current for the kink state decreases with 
I712I as shown in Fig.[5Jd), since the kink state gradually loses 
its stability when ly^l increases as discussed above. At the 
critical current, we observe a phase slip with a voltage pulse, 
and finally the system reaches a uniform state. 

The phase kink in Eq. |4]) between two TRSB pair states of 
a three-component superconductor is stable because the sys- 
tem takes different states in the left and right domains, and one 
cannot transform one state to the other by adjusting supercon- 
ductivity at the domain wall. This kink is thus protected by 
symmetry and is very different from those with bulk TRS as 
in Eq. where the states in the left and right domains are es- 
sentially the same except for a common phase factor, as shown 
in Fig. [T](middle), and the system evolves into a uniform state 
by rotating the phase of domains globally at the instability of 
the kinks. 

Consequences of TRSB - We have shown that spontaneous 
magnetic flux and voltage pulse appear during the nonequilib- 
rium evolution of superconductivity phase when a phase kink 
becomes unstable. Here we discuss possible experimental ob- 
servations on stable phase kinks at equilibrium. In the pres- 
ence of phase kink, the TRS is violated at the domain wall, 
namely <p\ - <p2 + or n. The variation of the superconductiv- 
ity amplitude is coupled with that of superconductivity phase, 
which can be checked by expanding the interband coupling 
term yyA,-A j cos(0; - <f>j) to the quadratic order of phase dif- 
ference. When the superconductivity is suppressed locally by 
nonmagnetic impurities, proximity effect at sample edge or 
heating, variations in superconductivity phases are induced, 
which in turn excites supercurrent and magnetic flux, as con- 
firmed numerically. 
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FIG. 3. (color online), (a) Suppression of the amplitude of super- 
conductivity at the domain wall when temperature denoted by a is 
increased in a two-band superconductor with y 12 = -0.9. (b) Phase 
diagram for the stability of phase kink, (c) Structure of the phase 
kink in the presence of supercurrent. Here aj = -7, y 12 = -0.5 
and the supercurrent J s = 2.94. (d) Stability of the phase kink upon 
current injection, with aj = -7, fij = 1 and nij = 2. 



FIG. 4. (color online). Numerical results of the magnetic field dis- 
tribution: (a) a two-band superconducting strip with phase kink in 
contact with a normal metal; (b) a three-band superconductor with 
TRSB in contact with a normal metal; (c) a three-band superconduc- 
tor with TRSB with an impurity. For(a),a, = -20, 712 = = 1. 
mi = 1 and ni2 = 3, with proximity lengths pjj = 2 and p 7 -# = 00; for 
(b), a, = 0,/3j = mj = 1, y 12 = 1, 713 = 1.2 and y 2 3 = 1-5, p n = 1 
and pj^i = 00; for (c), the same as (b) except for aj = 0.5 inside the 
impurity area, and p jt = 00. 



We study the proximity effect between a superconducting 
strip and a normal metal when a phase kink is present in the 
superconductor. In order to describe the proximity effect cor- 
rectly, a boundary condition between a multi-band supercon- 
ductor and a normal metals should be formulated. The bound- 
ary condition in terms of the Usadel equation has been derived 
in Ref. \29\. In the framework of phenomenological GL the- 
ory, the boundary condition for a single-band superconductor 
can be generalized straightforwardly to a multi-band one|22| 




where the off-diagonal coefficient p* with j + k accounts for 
the interband coupling while the diagonal coefficient j = k 
represents suppression of superconductivity as a consequence 
of the leakage of Cooper pairs at the interface. We minimize 
the GL energy numerically, and the results are presented in 
Fig. [4] (a). Spontaneous magnetic field is induced at the in- 
terface between the normal metal and superconductor at the 

gisition of the domain wall, which is strong enough (in Fig. 
a), H ~ \Qr 5 H C 2) to be measured experimentally by scan- 
ning SQUID, Hall, or magnetic force microscopy. The mag- 
netic field has opposite signs at the two interfaces, leaving a 
zero integration over the sample. 

In TRSB superconductors, stable domain walls associated 
with the variation of superconductivity phase and amplitude 
can be created by local perturbations, because the phase 
is coupled with the amplitude when TRS is violated. In 
Fig. |4|b), we consider the proximity effect between a three- 
band superconductor with TRSB and a normal metal. Mag- 
netic fluxes appear at the corners of the superconductor, as- 



sociated with sharp changes of phase gradients. In Fig. |4jc), 
we introduce an impurity by modifying a, locally. We see that 
magnetic flux is induced around the impurity. For supercon- 
ductors with TRS, no magnetic filed can be induced by the 
proximity effect or impurities, which implies a possible way 
to detect the TRSB in experiment. 

Discussions - In multi-band superconductors with Josephson- 
like interband coupling, the phase kinks as topological excita- 
tions can exist because of the multiple degenerate energy min- 
ima associated with the interband Josephson coupling. The 
phase kink suppresses the superconductivity nearby depend- 
ing on the ratio of the width of the phase kink to the super- 
conducting coherence length, which causes instability when 
the suppressed superconductivity is insufficient to maintain 
the phase coherence for the phase kink. The existence of the 
phase kink does not require breaking of additional symmetry 
besides t/(l). Instead, the presence of phase kink violates the 
TRS locally. When bulk multi-band superconductors break 
TRS, a new type of phase kink can be formed between the 
two TRSB pair states. The topological solutions (phase kinks) 
in multi-band superconductors discussed in the present work 
are different from the topological superconductors realized in 
materials with strong spin-orbit couplings. 

In ID, both phase kink in a superconductor with TRS and 
that between TRSB pair states are stable. In 2D, because of 
the attraction between opposite kinks, the domain wall col- 
lapses and the system reaches a uniform state. This is in 
accordance with the Derrick's theorem[30|, i.e. for an infi- 
nite system, the kink state is only stable in ID. Kinks can 
be pinned by the pinning centers where the superfluid densi- 
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ties are small, since the loss of superconductivity condensa- 
tion energy can be reduced by adapting the domain wall to 
the pinning centers. This may prevent the domain wall from 
collapsing and stabilize the kink in 2D. The kink can also be 
stabilized when vortices are present as discussed by Garaud et 
q/. |[T9l . Domain walls created by local heating or impurities 
in superconductors with TRSB are stable in 2D and 3D since 
they are enforced by external perturbations. 

Let us discuss the realization of the phase kink in Eq. 
In iron pnictide superconductors, interband scatterings 
are strong! 3 1|, an d thus the phase kinks are unlikely real- 
ized. While for the well-known two-band superconductor 
MgB2 and V^Si, it is revealed that interband scatterings are 
weak lF32l [33], which may allow for the excitation of stable 
phase kinks in low-temperature region. While for the realiza- 
tion of kink in Eq. Q, we need multi-band superconductors 
with TRSB. As discussed in Ref. [18|, the TRSB state can 
be achieved by chemical doping in iron pnictide supercon- 
ductors. Moreover, phase kinks can also be realized in hy- 
brid structures with two superconducting films coupled with 
Josephson coupling. Recently, Vakaryuk et al. proposed to 
realize the phase kinks in superconductors with the s± pairing 
symmetry by exploiting the proximity effect to a conventional 
s- wave superconductor. l34l 
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